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A  State  Space  Approach  to  Robustness  Analysis  and 
Synthesis  for  Nonlinear  Uncertain  Systems 

Wei-Min  Lu*  and  John  C.  Doyle'*' 


Abstract 

A  state-space  characterization  of  stability  and  performance  robustness  analysis  and  synthesis 
with  some  computationally  attractive  properties  for  nonlinear  uncertain  systems  is  proposed. 
The  robust  stability  and  robust  performances  for  a  class  of  nonlinear  systems  subject  to  bounded 
structured  uncertainties  are  characterized  in  terms  of  various  types  of  nonlinear  matrix  in¬ 
equalities  (NLMIs),  which  are  natural  generalizations  of  the  linear  matrix  inequalities 
(LMIs)  that  appear  in  linear  robustness  analysis.  As  in  the  linear  case,  scalings  or  multipliers 
are  used  to  find  storage  functions  that  give  sufficient  conditions  for  robust  performances;  these 
are  also  necessary  under  certain  assumptions  about  smoothness  of  the  storage  functions  and 
structure  of  the  uncertainty.  The  resulting  NLMIs  yield  convex  optimization  problems.  Unlike 
the  linear  case,  these  convex  problems  are  not  finite  dimensional,  so  their  computational  bene¬ 
fits  are  far  less  immediate.  Sufficient  conditions  for  the  solvability  of  robust  synthesis  problems 
are  developed  in  terms  of  NLMIs  as  well.  Some  aspects  of  the  computational  issues  are  also 
discussed. 


1  Introduction 

In  this  paper,  a  state-space  characterization  of  robust  stability  and/or  robust  performances  for  a 
class  of  nonlinear  systems  subject  to  bounded  structured  dynamic  uncertainties  is  proposed;  both 
analysis  and  synthesis  problems  are  addressed.  The  basic  block  diagram  for  an  uncertain  system 
is  as  follows, 


A 

G 

u 


w 


where  A  is  the  uncertainty  which  is  represented  as  a  nonlinear  time  invariant/varying  causal  op¬ 
erator  with  £2-gain  bounded  by  1,  G  is  the  nominal  system  which  is  nonlinear  time-invariant,  w  is 
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some  external  disturbance  vector,  and  z  is  the  regulated  signal  vector.  It  is  assumed  that  the  inter¬ 
connection  for  the  uncertain  system  is  well-posed  for  each  admissible  uncertainty.  The  robustness 
analysis  is  to  determine  that  under  what  conditions  for  nominal  system  G,  the  uncertain  system 
is  stable  and/or  satisfies  some  performance  for  all  admissible  uncertainty  A;  while  the  robustness 
synthesis  problem  is  to  decide  under  what  conditions  there  are  feedback  control  laws  for  the  uncer¬ 
tain  systems  such  that  the  closed  loop  uncertain  systems  have  the  required  robustness,  and  then 
design  the  control  law. 

The  systematic  treatment  of  robustness  analysis  for  such  uncertain  systems  can  be  traced  back 
to  at  least  the  60’s;  A  general  sufficient  condition  for  robust  stability  analysis  is  the  small  gain 
condition  where  the  £2- gain  of  the  nominal  system  G  is  bounded  by  1,  i.e.,  ||Gr||£2  <  1  [25,  36] . 
In  the  case  where  the  dynamic  uncertainty  is  unstructured,  i.e.  it  consists  just  one  full  block,  if 
both  the  uncertainty  and  the  nominal  system  are  linear;  or  the  uncertainty  is  linear  time-varying 
or  nonlinear,  and  the  nominal  plant  is  nonlinear  with  fading  memory,  then  the  small  gain  condition 
is  also  necessary  for  robust  stability  [10,  7,  27,  30],  where  the  necessity  means  in  case  ||G||£2  >  1, 
there  exists  an  admissible  destabilizing  perturbation  A. 

When  the  uncertainty  A  is  structured,  i.e.  it  consists  of  multiple  uncertainty  blocks,  A  = 
block-diag{A1?  A2,  •  •  ■,  Ajy},  the  sufficient  small  gain  condition  can  be  arbitrarily  conservative  for 
robust  stability.  In  the  case  where  both  uncertainty  and  nominal  system  are  linear  time  invariant, 
a  necessary  and  sufficient  condition  for  robust  stability  is  the  nominal  system  has  small  “structured 
gain”,  i.e.  its  structured  singular  value  //A(G)  <  1  [9].  A  conservative,  sufficient  condition  for  robust 
stability  [9,  4]  in  this  case  is  that  the  nominal  system  has  scaled  small  gain:  || DG'D~1||Z.2  <  1  for 
some  scale  D  commuting  with  A.  Under  certain  assumptions,  this  condition  is  also  necessary, 
such  as  when  the  the  nominal  plant  G  is  linear  time  invariant  and  the  uncertainty  is  allowed 
to  be  time- varying  [28,  14,  8].  Furthermore,  if  G  is  finite  dimensional,  the  computation  of  the 
condition  \\DGD~1\\C^  <  1  can  be  converted  to  a  finite  dimensional  linear  matrix  inequality 
(LMI),  which  is  computationally  attractive,  (see  [6]  for  a  tutorial  review  of  LMIs  and  their  use.) 
In  the  case  where  the  nominal  system  is  nonlinear  with  fading  memory,  the  uncertainty  structure 
is  linear  time- varying  or  nonlinear,  the  scaled  small  gain  is  necessary  as  well  [26],  although  the 
computational  aspects  are  not  as  well  resolved. 

As  for  robust  performance  analysis  problem,  since  the  performance  robustness  analysis  problem 
can  be  transferred  into  a  robust  stability  problem  with  structured  uncertainty  by  adding  an  extra 
“uncertainty”  block  [11],  the  above  small  gain  arguments  still  apply  in  this  case;  see  [11,  21,  24, 14,  8] 
and  references  therein  for  this  consideration. 

The  above  analysis  of  robust  stability  and  performance  for  uncertain  systems,  which  are  treated 
in  the  input/output  setting,  are  essentially  reduced  to  (scaled)  system  gain  analysis.  Therefore, 
it  is  possible  that  the  internal  behaviors  of  the  systems  are  ignored  during  the  analysis  in  the 
input/output  setting.  In  addition,  the  computational  implications  for  those  characterizations  are 
not  clear  for  the  nonlinear  case.  These  concerns  can  be  remedied  by  taking  a  state  space  treatment. 
A  basic  internal  consideration  in  the  state  space  analysis  is  to  see  if  the  system  is  asymptotically 
stable.  For  linear  systems,  the  asymptotic  stability  is  guaranteed  by  the  input/output  stability 
if  the  system  is  detectable  and  stabilizable;  and  there  are  systematic  ways  to  get  such  internally 
well-behaved  realization.  In  particular,  the  £2~stability  and  -performance  robustness  for  linear 
systems  can  be  easily  characterized  in  state-space  in  terms  of  the  results  of  -performance  analysis 
in  state  space  [12,  19,  22,  23,  16,  21,  1,  6].  However,  the  story  for  state-space  treatment  of  the 
nonlinear  systems  is  different,  since  the  small  gain  theorem  can  only  guarantee  the  system  to 
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have  bounded-input/bounded-output  property;  it  does  not  reveal  internal  information  about  the 
systems,  and  there  is  more  required  to  insure  asymptotic  stability.  Moreover,  the  performance  and 
robust  stability  problems  are  essentially  different,  although  the  characterizations  are  similar. 

In  this  paper,  we  give  state-space  characterizations  of  stability  and  performance  robustness  for 
nonlinear  uncertain  system,  and  consider  both  analysis  and  synthesis  problems.  By  robust  stability, 
we  mean  that  the  feedback  system  is  asymptotically  stable  for  each  admissible  uncertainty;  the 
robust  performance  means  that  the  uncertain  system  is  asymptotically  stable  and  has  C2- gain  <  1. 
The  treatments  of  the  robustness  issues  in  this  paper  are  motivated  by  the  small  gain  theorem 
and  its  recent  extensions,  together  with  the  LMI  characterization  of  results  in  the  linear  case. 
Essentially,  one  of  the  main  results  in  this  paper  implies  that  under  some  additional  (stabilizing) 
conditions,  if  both  the  (structured)  uncertainty  and  the  (scaled)  nominal  system  are  bounded  by 
1  in  £2  sense,  then  the  uncertain  system  is  robustly  stable  or  has  robust  performance  and  we 
characterize  all  of  the  conditions  in  terms  of  nonlinear  matrix  inequalities  (NLMIs). 

Although  other  characterizations,  such  as  Hamilton -Jacobi  equations/inequalities,  exist  and  are 
equivalent  to  the  NLMIs  in  this  paper,  the  characterizations  in  terms  of  NLMIs  offer  potentially 
attractive  computational  features.  In  particular,  like  the  linear  case,  the  NLMIs  trivially  give  convex 
conditions  on  the  unknowns.  Unfortunately,  the  NLMI  conditions  involve  neither  a  finite  number  of 
unknowns  nor  a  finite  number  of  constraints,  so  the  computational  advantages  are  far  less  immediate 
than  for  LMIs.  Clearly  much  additional  work  will  be  needed  on  the  computational  aspects  and 
sophisticated  approximation  techniques  may  be  required  to  make  the  NLMI  computation  feasible. 

The  rest  of  the  paper  is  organized  as  follows.  In  section  2,  some  results  for  asymptotic  stability 
and  £2-gain  analysis  for  nonlinear  systems  are  reviewed,  and  these  characterizations  are  reformu¬ 
lated  as  NLMIs.  In  section  3,  the  stability  robustness  analysis  is  conducted;  both  unstructured  and 
structured  uncertainty  cases  are  treated  independently,  although  the  unstructured  uncertainty  is  a 
special  case  of  the  structured  uncertainty;  the  former  is  used  to  motivate  the  latter.  In  section  4, 
the  robust  performance  analysis  is  conducted;  the  general  structured  uncertainty  case  is  considered, 
and  two  cases  where  the  plants  are  generally  causal  and  strictly  causal  are  dealt  with  respectively. 
The  characterizations  of  both  robust  stability  and  robust  performance  are  in  terms  of  NLMIs.  In 
section  5,  we  deal  with  the  robustness  synthesis  problem;  we  just  take  the  state  feedback  perfor¬ 
mance  robustness  synthesis  problem  as  an  example,  and  the  solvability  conditions  are  also  in  terms 
of  NLMIs.  We  address  some  computational  issues  for  robustness  analysis  and  synthesis  in  section 
6.  In  the  concluding  remarks,  we  briefly  discuss  the  scaled  small  gain  theorem  for  bounded-input- 
bounded  output  stability  of  nonlinear  uncertain  systems,  which  is  the  motivation  for  the  treatments 
in  this  paper. 

Conventions 

The  following  conventions  are  made  in  this  paper.  R  is  the  set  of  real  numbers,  R+  :=  [0,  00)  C  R. 
R”  is  n-dimensional  real  Euclidean  space;  ||-||  stands  for  the  Euclidean  norm.  X  (or  X,:)  is  the 
state  set  which  is  a  convex  open  subset  of  some  Euclidean  space  and  contains  the  origin.  Rnxm 
((Qrcxm)  ^jle  ge^  an  n  X  m  real  (complex)  matrices.  The  transpose  of  some  matrix  M  £  Mnxn 
is  denoted  by  MT .  By  P  >  0  (P  >  0)  for  some  Hermitian  matrix  P  £  Rnxn  or  (Cixm)  we  mean 
that  the  matrix  is  (semi-)positive  definite.  A  function  is  said  to  be  of  class  if  it  is  continuously 
differentiable  k  times;  so  C°  stands  for  the  class  of  continuous  functions.  A  function  V  :  X— ^R+  is 
positive-definite  if  V(x )  >  0,  V(x )  =  0=>x  =  0,  and  limipn^oo  V(x )  =  00  on  X.  A  C°  matrix- valued 
function  P  :  X— >R"X"  is  positive  definite  if  P(x)  is  positive  definite  for  each  a;  £  X  and  there  exist 
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a  >  0  such  that  P(x)  >  al  for  all  x  £  X.  £2(R+)  stands  for  the  function  space  which  is  defined 
as  the  set  of  all  vector-valued  functions  u(t)  on  R+  such  that  ||rz||2  :=  (Jo°°(||w(f)||2  dt)1!2  <  oo,  and 
£^(R+)  is  its  extended  space  which  is  defined  as  the  set  of  the  vector- valued  functions  u(t)  on  R+ 
such  that  PTu(t )  £  £2(R+)  for  all  T  £  R+,  where  PT  is  the  truncation  operator. 


2  Preliminaries 

In  this  section,  we  will  review  some  standard  results  about  stability  and  'Hoo -performances  of  a 
class  of  nonlinear  time  invariant  systems.  We  will  reformulate  these  characterizations  in  terms  of 
so-called  nonlinear  matrix  inequalities  (NLMIs). 

The  system  considered  in  this  section  has  the  following  control  affine  realization, 


G  : 


%  =  f{x)  +  g(x)w 
z  =  h(x )  +  k(x)w 


(1) 


where  x  G  Rra  is  state  vector,  w  G  Rp  and  z  £  l!  are  input  and  output  vectors,  respectively.  It 
is  assumed  that  f,g,h,k  G  C°  are  verctor  or  matrix  valued  function,  and  /( 0)  =  0,h(0)  =  0. 
From  now  on  we  will  assume  the  system  evolves  on  a  convex  open  subset  X  C  Rn  containing  the 
origin.  Thus,  0  G  Rn  is  the  equilibrium  of  the  system  with  w  =  0.  The  state  transition  function 
4>  :  R+  x  X  X  £2(^+)  — ►  X  is  so  defined  that  x  =  <^(jT,  x0,  w*)  means  that  system  G  is  driven  from 
initial  state  x0  to  state  x  in  time  T  by  the  control  action  w*  G  £2(R+). 

Note  that  in  many  cases  system  (1)  can  be  rewritten  (nonuniquely)  as  the  following  form  which 
is  also  used  in  this  paper. 


G  : 


x  =  A(x)  x  +  B(x)w 
z  =  C(x)x  +  D(x)w 


(2) 


where  x  G  Mn  is  state  vector,  w  G  Rp  and  z  £  K9  are  input  and  output  vectors,  respectively.  We 
will  assume  A,  B,  C,  D  are  C°  matrix-valued  functions  of  suitable  dimensions. 

Definition  2.1  (i)  The  system  (1)  or  (2)  is  reachable  from  0  if  for  all  x  G  X,  there  exist  T  G  R+ 
and  w*(t )  G  £2[0,T]  such  that  x  =  <f>(T,0,w*); 

(ii)  The  system  (1)  or  (2)  is  (zero-state)  detectable  if  w  =  0  and  for  all  x  G  X,  z(t)(= 
h(<f){t,  x,  0)))  =  0  =>  4>(t,  a:,0)-^0  as  t— »oo. 

2.1  Asymptotic  Stability 

Definition  2.2  Consider  the  nonlinear  dynamical  system  G.  It  is  asymptotically  stable  if 
lim  (p(t,  x0 , 0)  =  0 

t-+  OO 

for  any  initial  state  x0  G  X. 

We  have  the  following  results  about  asymptotic  stability. 
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Proposition  2.3  (i)  (Lyapunov)  The  equilibrium  0  of  the  system  G  (1)  with  w  =  0  is  asymptoti¬ 
cally  stable  if  and  only  if  there  exists  a  C1  positive  definite  function  V  :  X— such  that  V(x)  <  0, 
where  V(x)  =  0  if  and  only  if  x  —  0. 

(ii)  (LaSalle)  The  equilibrium  0  of  the  system  G  (1)  with  w  —  0  is  asymptotically  stable  if  and 
only  if  there  exists  a  C1  positive  definite  function  V  :  X— >-M+  such  that  V{x)  <  0,  and  those  x(t) 
with  x  =  f(x)  such  that  =  0  satisfy  a:(t)— >-0  as  t^oo. 

The  positive  definite  functions  V  in  the  above  theorem  is  called  Lyapunov  functions.  Note 
that  the  above  Lyapunov  theorem  can  be  explicitly  restated  as  follows. 

Theorem  2.4  Consider  system  (1).  It  is  asymptotically  stable  around  0,  if 

(i)  there  is  a  C1  positive  definite  function  V  :  X— ^R+  such  that 

x)f(x )  <  0.  (3) 

for  all  x  €  X  \  {0}  with  x  ^  0;  or 

(ii)  (1)  is  detectable,  and  there  is  a  C1  positive  definite  function  V  :  X— >R+  such  that 

j^(x)f(x)  + hT(x)h(x)  <  0.  (4) 

for  all  x  e  X. 

Remark  2.5  It  is  noted  that  the  conditions  of  Theorem  2.4  for  stability  are  convex  conditions. 
That  is,  in  both  cases,  the  positive  solutions  V  satisfying  (3)  and  (4)  form  convex  sets,  respectively. 
This  trivial  fact  has  only  been  exploited  systematically  in  the  linear  case,  but  we  hope  that  numerical 
techniques  may  be  developed  to  exploit  it  in  the  nonlinear  case  as  well. 


2.2  £2-Gains  of  Nonlinear  Systems 

Definition  2.6  The  system  G  (1)  with  initial  state  x(O)  =  0  is  said  to  have  £2-gain  less  than 
or  equal  to  7  for  some  7  >  0  if 

[  IKt)||2dt  <  72  f  IK<)ll2rft  (5) 

Jo  Jo 

for  all  T  >  0  and  w(t )  £  £2(M+),  and  z(t )  =  h{x(t))  +  k{x(t))w(t )  with  x(t)  —  <f>(t,0,w(t)). 

The  following  results  characterizes  £2-gains  for  a  class  of  nonlinear  systems  which  are  asymptoti¬ 
cally  stable  in  terms  of  NLMIs.  The  reader  is  referred  to  [18]  for  more  about  NLMI  characterizations 
of  £2-gains  for  nonlinear  systems. 


Theorem  2.7  Consider  system  G  given  by  (1))  with  R(x )  =  I—kT(x)k(x )  >  0,  it  is  asymptotically 
stable  and  has  C2-gain  <  1  if 

(i)  [f(x),h(x)]  is  detectable,  and  there  exist  a  C1  positive  definite  function  V  :  X— >R+  such 

that 


f 7(x)f(x)  +  hT{x)h(x)  \%{x)  g(x)  +  hT  (x)k(x) 

2 9T(x)~^(x)  +  kT(x)h(x)  kT(x)k{x)  -  I 
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(6) 


for  all  x  £  X. 

(ii)  there  exist  a  C1  positive  definite  function  V  :  X— such  that 

| j(x)f(x)  +  hT(x)h(x)  || j(x)g(x)  +  hT(x)k(x)  1 

_|  gT(x)^f(x)  + kT(x)h(x)  kT(x)k(x)-I  J 

for  all  x  £  X  \  {0}. 


Proof  By  Schur’s  complement  argument,  we  have  that  (6)  is  equivalent  to 
‘ H(V,x )  :=  ~(x)f(x)  +  hT(x)h(x)+ 

+  hT {x)k(x))(I -  kT(x)k(x)y1(x)(^gT(x)^^(x)  +  kT(x)h(x))  <  0.  (8) 

which  is  a  Hamilton-Jacobi  inequality  [31],  and  implies 


dV_ 

dx 


(x)f(x)  <  -hT(x)h(x)+ 

+  hT (x)k(xW ~  kT(x)Hx))~\x)(\gT(x)^r(x)  +  kT{x)h{x)). 


Thus, 


dV . 


V(X)  =  -gy(x)(f(x)  +  9(x)w) 

dv, 


dx 

1  BV  nn  1  BVT 

•(-2  ~dx^9^  +  h  -  kT(x)k(x)y1(x)(2^T^^y^  +  kT(x)Kx)) 


( x)g(x)w  -  hT(x)h(x)+ 


<  IKOH2  -  ||*(<)||a  - 


1  f)VT 

R1/2(x)w(t)  -  R~1/2(x)kT(x)h(x)  -  -R~1/2(x)gT (x)-^—(x) 

£  (J  Jb 


The  latter  inequality  follows  by  completion  of  squares  for  (2.2).  Therefore, 

*"(*)- (IKOII2  -  ||*(<)ll2)  <  0. 

Take  the  integral  from  t  =  0  to  t  =  T,  the  above  inequality  implies  that  the  system  has  £2-gain 
<  1  since  V(x)  >  0. 

On  the  other  hand,  note  that  H(V,  x)  <  0  implies 
dV 

-fa  (x)f(x)  +  hT(x)h(x)  <  0; 

and  [f(x),h(x)]  is  detectable  by  the  assumption,  then  the  stability  is  confirmed  by  theorem  2.4. 

(ii)  By  using  Schur  complement  argument,  we  have  that  (7)  is  equivalent  to  H(V,x)  <  0  for  all 
x  G  X  \  {0},  where  H  is  defined  in  (8).  It  follows  from  the  similar  arguments  as  above  that 

*(*)<  IK*)lla-lkWlla  +  w(v,*)  (9) 
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<  IH*)II2  -  \W) 


It  is  easy  to  see  that  the  latter  inequality  implies  the  £2-gain  <  1.  Now  take  w(t)  =  0,  (9) 
becomes 


V(x)  <  -  \\z(t)\\2  +  H(V,x) 


So  V(x)  =  0  implies  H(V,x)  =  0,  which  in  turn  implies  x  =  0. 
Lyapunov  function,  and  system  G  is  asymptotically  stable. 


Therefore,  V  :  X^R+  is  a 

□ 


Remark  2.8  Note  that  conditions  in  Theorem  2.7  are  affine  in  V(»),  and  all  such  solutions  form 
convex  sets.  These  inequalities  are  actually  differential  linear  (or  affine)  matrix  inequalities,  but  we 
will  refer  to  them  as  nonlinear  matrix  inequalities  (NLMIs)  to  emphasize  their  use  in  nonlinear 
problems.  All  of  the  conditions  that  are  derived  for  the  analysis  problems  in  the  remainder  of  this 
paper  are  similarly  convex,  and  this  property  will  not  be  discussed  for  each  problem. 


2.3  Further  Remarks:  Alternative  Characterizations 


It  may  be  helpful  for  computation,  especially  in  the  synthesis  problem,  to  have  some  alternative 
characterizations,  instead  of  using  (3),  (6),  and  (7);  however,  the  price  will  often  be  paid  for 
increasing  conservatism.  For  example,  as  in  [18],  system  (2)  is  considered  and  (7)  is  replaced  by 
the  following  NLMI  in  the  characterization, 

'  AT (x)P(x)  +  P{x)A(x)  +  CT(x)C{x )  P{x)B{x)  +  CT(x)D(x) 

Bt(x)P(x)  +  Dt{x)C(x)  Dt(x)D(x)-I 


<  0.  (10) 


for  some  positive  definite  C°  matrix  valued  function  P  :  X^RnXn  such  that  §7(2)  =  2xTP(x)  for 
some  C1  function  on  X  with  V(0)  =  0.  In  fact  if  such  function  P  exists,  (10)  implies 


xT  0 

AT(x)P(x )  +  P{x)A(x)  +  CT(x)C(x ) 

P(x)B(x)  +  CT(x)D(x)  ' 

x  0 

0  I 

BT(x)P(x)  +  DT{x)C(x) 

Dt{x)D(x)-I 

0  I 

for  all  x  6  X  \  {0},  which  is  exactly  (7).  We  will  see  the  characterizations  in  terms  of  (10)  have 
some  computationally  attractive  properties. 

It  is  noted  that  the  NLMI  (10)  has  positive  definite  solution  P  :  X->RnX”  implies  that  system 
(2)  point-wise  has  £2- gain  <  1;  however,  to  guarantee  the  performance  for  the  original  system, 
it  is  additionally  required  that  there  exists  a  function  V  :  X-^R  such  that  j^-(x)  =  2 xTPT(x). 
Lemma  6.5  provides  a  characterization  of  a  class  of  matrix-valued  functions  P  :  X-^Rnxn  which 
satisfy  this  additional  requirement.  The  same  arguments  are  also  true  for  the  stability  analysis. 

To  see  the  conservativeness  of  NLMI,  such  as  (10),  it  may  be  noted  that  even  its  point-wise 
solution  is  not  necessary  for  either  stability  or  performance,  as  is  well-known.  This  is  shown  in  the 
following  example. 


Example  2.9  Consider  the  following  system  of  order  two, 

{*1  =  —X\  +  xxx\ 

x2  =  -x\x2  -  x2 
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which  evolves  on  R2.  Take  a  positive  definite  (quadratic)  function  V{x1,x2)  =  x\  +  x\,  which  is 
actually  a  Lyapunov  function  for  the  system  and  satisfies  (3).  In  fact,  V(x1,x2)  =  —  2(x2  +  x\), 
which  is  negative  definite;  the  system  is  therefore  asymptotically  stable. 

On  the  other  hand,  the  system  can  be  rewritten  as  the  form  (2),  i.e.,  x  =  A( x)x  for  x  E  R2. 
There  are  two  representations  in  which 


A(x) 


1  +  ^2  0 
0  -l-x\ 


and  A(x ) 


-1  xxx2 
-xxx2  -1 


Obviously,  neither  A(x)  ’s  are  point-wise  asymptotically  stable  on  R2.  This  may  be  compared  to  the 
pointwise  use  of  (10),  which  is  discussed  further  in  Section  6.3. 


3  Stability  Robustness  of  Uncertain  Systems 

Consider  the  following  uncertain  system  which  is  described  as  a  feedback  system  set, 


y 


u 


where  the  nominal  system  G  has  the  realization  similar  to  (2),  i.e., 
G  : 


x  =  f{x)  +  g{x)u 
y  =  h(x)x  +  k(x)u 


and  the  uncertainty  A  is  fed  back  to  the  nominal  system  and  belongs  to  a  set: 

A  :=  {A  :  A  is  a  causal  asymptotically  stable  nonlinear  time-varying  system} 


(11) 


(12) 


In  this  section,  the  case  where  the  uncertainty  is  norm-bounded  is  considered.  In  particular,  it 
is  assumed  that  A  is  in  the  following  uncertainty  subset  of  A. 


BA  :=  {A  :  A  G  A  and  has  C2- gain  <  1} 


(13) 


Definition  3.1  The  uncertain  system  is  robustly  stable  if  for  each  A  E  BA,  the  feedback  system 
is  well-posed  and  asymptotically  stable  around  0. 


From  now  on,  it  is  always  assumed  that  each  uncertainty  A  E  A  is  fed  back  to  the  nominal 
system  G  in  the  well-posed  manner.  In  this  section,  we  examine  the  robust  stability  in  two  cases 
where  the  uncertainties  A  are  unstructured  and  structured,  respectively. 
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3.1  Systems  with  Unstructured  Perturbations 

We  first  examine  the  uncertain  structure  A  in  a  simple  case,  i.e.  A  has  the  following  input  affine 
realization. 

A  .  f  £  =  fd(0  +  9d(0v 
\u  =  MO 

with  fd( 0)  =  0  and  hd( 0)  =  0.  It  is  known  that  A  has  £2-gain  <  1  if  and  only  if  there  exists  a  well 
defined  storage  function  U  :  X0— R+  for  this  system  with  respect  to  the  supply  rate  ||y||“  -  ||u||  (see 
[35]);  if  U(0  is  differentiable,  then  A  has  £2-gain  <  1  if  and  only  if  the  following  Ilamilton-Jacobi 
inequality  holds  (see  [31]): 

m  t=  ~(OUO  +  )?£(()  + <  o 

Then  it  can  be  shown  that  the  above  inequality  is  equivalent  to 

*7(0  <  II»IIs-INI2  +  ^(0 

for  all  possible  y(t).  Motivated  by  this  observation,  we  make  the  following  assumption  for  a  general 
uncertain  nonlinear  structure  A  £  BA. 


Assumption  3.2  For  each  A  £  BA,  it  has  the  following  realization: 

f  £  =  fd(t,t,v) 

\  u  =  hd(£,t,y) 

which  evolves  on  X0  and  £  =  /d(£,  t,  0)  is  asymptotically  stable  at  0  £  X0;  in  addition,  there  exists 
a  C1  storage  function  U  such  that  U{()  <  J  -  ||u||2  +  with  some  negative  definite  function 
$  :  X;;-E+. 

We  have  the  following  theorem  regarding  the  robust  stability. 


Theorem  3.3  Under  assumption  3.2,  the  uncertain  system  is  robustly  stable  if  there  is  a  positive 
definite  C1  function  V  :  X— »R+  satisfying  the  following  inequality: 

%(x)f{x)  +  hr(x)h(x)  \{x)^g{x)  +  hT  (x)k(x)  1  (  , 

.  \9T(x)^Ur(x)  +  kT(x)h(x)  kT(x)k{x)-I  J 

for  all  x  £  X\{0}. 


Proof  Consider  (14),  note  that  that  R(x)  :=  I  —  kT (x)k(x)  >  0,  by  using  Schur  complement 
argument,  we  have  that  it  is  equivalent  to  the  following  Hamilton-Jacobi  inequality: 

f)V 

H(V,  x )  :=  —  (x)f(x)  +  hT(x)h(x)+ 


+(^J^(x)9(x)+hT(x)k(x))(I-kT(x)k(x))  1(x)(^gT(x)^-(x)+kT(x)h(x))  <  0.(15) 
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for  all  x  G  X\{0}.  Take  V  as  defined  in  the  statement,  similar  argument  to  the  proof  of  Theorem  2.7 
yields  that 

V(x)<\\u(t)\\2-\\y(t)\\2  +  n(V,x).  (16) 

On  the  other  hand,  by  assumption  3.2,  for  each  A  E  BA,  there  is  a  positive  definite  function 
U  :  X0-^R+,  such  that 

u( o  <  ||y(<)ll2  -  IK*)II2  +  VKO-  (17) 

for  some  negative  definite  function  if  on  X0,  where  £  is  the  state  vector  of  A  defined  on  X0. 

Next,  define  a  positive  definite  function  W  on  X  x  X0  as 

W(x,0  =  V(x)  +  U(0  (18) 

So  from  (16)  and  (17),  it  follows  that 

W(x,£)<H(V,x)  +  iP(O<0 

Thence,  if  W(x,(f)  =  0,  then  H{V,x)  =  0  and  if(£)  =  0;  the  former  condition  implies  x  =  0 
by  (15),  and  the  latter  implies  £  =  0  by  assumption  3.2.  Thus,  W  :  X  X  X0— >-R+  is  a  Lyapunov 
function  for  the  feedback  system,  and  the  system  is  asymptotically  stable.  Therefore,  the  uncertain 
system  is  robustly  stable. 

□ 


3.2  Systems  with  Structured  Perturbations 

In  this  subsection,  we  assume  the  uncertainty  has  the  following  structure: 

BA  :=  {A  =  block-diag{Ai,  A2,  •  •  •,  Ajv}  :  A  G  A  has  £2-gain  <  1}.  (19) 

Note  that  A  :=  block-diagjAi,  A2,  •  •  ■ ,  Ajv}  G  BA  if  and  only  if  block  A,-,  which  is  a  nonlinear 
time- varying  and  causal  system,  has  £2-gain  <1.  As  in  the  previous  subsection,  the  following 
assumption  is  made. 

Assumption  3.4  For  each  A  :=  block-diag{ Ai,  A2,  •  •  Ajy}  E  BA,  A,-  (i  £  {1,2,  •  ■  -,N})  has 
the  following  realization: 

f  ii  =  Mti,  t,  Vi) 

(  Wj  —  h  i  ((("f ,  t .  ) 

which  evolves  on  X,  and  has  a  unique  asymptotically  stable  equilibrium  at  0  E  X,  for  yi  =  0;  in 
addition,  there  is  a  C1  storage  function  C7f  such  that  f/,-(^)  <  \\yi\\2  —  1 1 1 1 2  +  ^«(£;)  with  some 
negative  definite  function  if>i  :  Xj— >-R+. 

To  reduce  the  conservatism  arising  from  the  uncertainty  structure,  we  can  perform  a  standard 
manipulation  as  in  the  following  diagram, 
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where  D  is  some  real  invertible  matrix.  Note  that  the  above  uncertain  system  is  the  same  as  the 
one  in  the  original  diagram.  Define  a  real  valued  matrix  set  V  as 

V  :=  {b\ock-dia,g{diI,d2I,- • -,dNI}  :  <£■  £  K, d,-  >  0}  (20) 

where  each  of  the  identity  matrices  is  compatible  with  the  corresponding  nonlinear  uncertainty  A,. 
It  is  noted  that  DA  =  AD  for  all  D  £  V  and  A  £  BA;  and  for  each  D  £  V,  A  £  BA  if  and  only 
if 

DAD~l  :=  {diA1df1,d2A2d21,-  •  -,dNANd]fi}  £  BA. 

In  fact,  the  £2-gains  of  A  and  DAD~l  are  the  same  if  D  £  V. 

Therefore,  we  can  view  DAD -1  as  a  legal  (transformed)  uncertainty  structure;  it  satisfies 
assumption  3.4  as  does  A.  Thus,  we  may  consider  the  scaled  system  DGD~l  which  reduces  the 
conservatism  arising  from  the  uncertainty  structure  (see  section  7  for  the  motivation).  We  have 
the  following  theorem  about  the  robust  stability  for  the  structured  uncertain  systems,  which  gives 
a  natural  NLMI  generalization  to  nonlinear  systems  of  the  LMI  conditions  for  the  linear  case. 


Theorem  3.5  Under  assumption  3-4,  the  structured  uncertain  system  with  nominal  plant  (11)  is 
robustly  stable  if  there  exist  a  positive  definite  C1  function  V  :  X— >R+  and  a  positive  definite  matrix 
Q  £  V  such  that  the  following  NLMI  holds: 


M(V,Q,x) 


| j(x)f(x)  +  hT(x)Qh(x)  §§£(*)$(*)  +  hT(x)Qk{x) 

\gT (x)^f(x)  +  kT(x)Qh(x )  kT (x)Qk(x)  -  Q 


(21) 


for  all  x  £  X  \  {0}. 


Proof  Consider  (21),  note  that 


I  0 

0  Q-1/2 


M(V,Q,x) 


I  0 
0  Q-1/2 


fjr(-T  )f(x)  +  hT(x)Qh(x) 

|  Q-1DgT(x)^f(x)  +  Q~l!2kT  (x)Qh(x)Q~1!2 


II 7(x)9(x)Q  1/2  +  hT(x)Qk(x)Q  1/2 
Q-1/2kT(x)Qk(x)Q^1D  -  I 


<  0. 


Define  g(x )  =  g(x)Q  1/l2,  h(x)  =  Q1/2h(x)  and  k(x)  =  Q1!2k(x)Q  1/2.  Using  Schur  comple¬ 
ments  argument,  We  have  that  the  above  inequality  is  equivalent  to  the  following  two  inequalities. 


R(x )  :=  I  -  kT(x)k(x )  >  0, 
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and 


H(V, Q,  x)  :=  —(x)f(x)  +  /iT(x)A(a;)+ 

H^^(x)g(x)+'hT(x)k(x))(I-kT(x)k(x))~\x)(^gT(x)^-(x)+kT(x)h(x))  <  0.(22) 
for  all  x  £  X  \  {0}.  Take  V  as  defined  in  the  statement,  and  define  u  —  Q1!2u  and  y  =  Q1!2y,  then 
V(x)  =  |jj(a0(/(a0  +  9{x)u) 


IK*)I 


\m\\2  - 


1  f)VT 

R1/2(x)u(t)  +  R~1/2(x)(kT(x)h(x)x  +  -^T(x)-^-(x)) 


+  'H  ( V,  Q ,  x ) . 


n*)<\m\\2  -\m\\2 +'h(v,q,x).  (23) 

On  the  other  hand,  notice  that  Q 1/2  E  V ;  denote  =  biock-diag{gi/,  q2I,  ■  ■  ■ ,  From 

the  assumption  3.2,  for  each  A  £  BA,  there  is  a  positive  definite  function  Ui  :  X  — >R+  for  nonlinear 
system  g,-  A,-?,"1  for  each  i  £  {1,2,--*,  iV}  such  that 

#<(&)  <  l|y;(*)H2  -  +  <M£«)  (24) 


for  some  negative  definite  function  on  X„  where  u,  =  <ji  U{ ,  f/,  =  g,-y,-  and  £,•  is  the  state  vector  of 
A;  on  X,-.  Note  that 


Ul 

’  Vi 

u2 

V2 

.  _ 

>  V  = 

.  y* . 

Therefore, 

imi2  =  ziii^ii2’  ii^ii2  =  o&n2 

i=l  i=l 

Next,  define  a  positive  definite  function  W  on  X  x  Xi  x  •  •  •  x  Xjv  as 

AT 

W{x,  Zu  ■  ■  • ,  Zn)  =  V{x)  +  Y,  Ui(Zi)  (25) 

2  =  1 

So  from  (23)  and  (24),  it  follows  that 

i v 

W^x,6,---,6v)  <  ||n(t)||2  -  P(t)||2  +  7f(F,g,x)  +  ^((|yi(t)||2  -  ||iii(t)||2  +  &(&)) 


N 

<H(P,x)  +  YMti) 

2  =  1 

Thence,  if  W(x,Zi,  -  •  •  ,Zn)  =  0,  then  H(V,Q,x)  =  0  and  i>i(Zi)  —  0  for  each  *;  the  former 
condition  implies  a:  =  0  by  (22),  and  the  latter  implies  &  =  0  by  the  assumption  3.4.  Thus, 
IT  :  X  X  Xi  X  •  •  •  x  Xjv-»R+  is  a  Lyapunov  function  for  the  feedback  system,  and  the  system  is 
asymptotically  stable.  Therefore,  the  uncertain  system  is  robustly  stable.  □ 
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3.3  More  about  Robust  Stability 

In  the  previous  discussion,  it  is  assumed  that  the  uncertain  system  is  well  posed  for  each  A  £  BA. 
It  is  known  that  the  well-posedness  is  guaranteed  when  any  one  of  A  and  G  is  strictly  causal 
[34],  In  this  section,  we  will  consider  the  case  when  G  is  strictly  causal.  Accordingly,  the  class  of 
uncertainty  BA  is  enlarged.  (Strictly  causal  uncertainty  case  can  be  considered  analogically). 

The  nominal  system  G  has  the  following  realization 
r  A  *  =  /(*)  +  s(x)u  font 

G  ■  \  y  =  />w  l2&) 

The  uncertainty  is  structured  as  follows 

BA  :=  {A  =  block-diag{A1,  A2,  ■  •  Ajv}  :  A  £  A  has  £2-gain  <  1}. 

The  following  assumption  is  made. 

Assumption  3.6  For  each  A  :=  {A  =  block-diag{ Ai,  A2,  •  •  ■ ,  AN}  £  BA,  A;  (i  £  {1, 2,  •  •  • ,  N} ) 
has  the  following  realization: 

^  'U'i  —  Vi) 

which  evolves  on  X,  and  &  —  fi(£i,t,  0)  is  asymptotically  stable  around  0  £  X,-;  in  addition,  there 
is  a  C1  storage  function  Ui  such  that  t/j(£;)  <  ||?/j||“  -  ||uj||2. 

Similarly,  we  consider  a  scaled  nominal  system  to  reduce  the  conservatism  which  arises  from 
the  structural  constraints  of  uncertainty.  The  scaling  matrix  set  V  is  defined  as  before 


V  \=  {block-diag{dx/,  d2I,  •  •  • ,  dNI}  :  d,  £  R,  d;  >  0} 


We  have  the  following  theorem  about  the  robust  stability  for  this  class  of  uncertain  systems 
where  the  nominal  systems  are  strictly  causal. 

Theorem  3.7  Under  assumption  3.6,  the  uncertain  system  with  nominal  system  as  in  (26)  is 
robustly  stable  if  there  exist  a  positive  definite  function  V  :  X— E+  and  a  positive  definite  matrix 
Q  £  V  such  that  the  following  NLMI  holds: 

| j{x)f(x)  +  hT(x)Qh(x) 

yT(x)^f{x) 

for  all  x  G  X  \  {0}. 


-Q 


<  o 


(27) 
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Proof  (27)  implies  that 

H(V,Q,x):=  ^(fc)/(s)  +  ^^(x)g(x)QgT(x)^-(x)  +  hT(x)Qh(x)  <  0 

for  all  x  £  X  \  {0}.  Take  V  :  X-4R+  and  U{  :  X{^R+  as  given.  Define  a  positive  definite  function 
W  on  X  x  Xi  x  •  •  •  x  Xjv  as 


N 

W(x,  &  •  •  • , &)  =  V(x)  +  E  (28) 

i=l 

Similar  arguments  used  in  the  section  3.2  yield  that 

Thence,  if  lT(x,6,  •  •  -,6v)  =  0,  then  H(V,Q,x )  =  0;  it  in  turn  implies  x  =  0. 

On  the  other  hand  x  =  0  implies  y  -  0,  so  the  feedback  system  evolves  on 

{(*,&,  ”-,60  6  X  x  Xi  x  •  •  •  x  Xiv  :  x  =  0,6  =  /<(&,*,  0)Vi  £  {1, 2,  •  •  •,  iV}}. 

But  by  the  assumption  3.6,  6  =  /»(&,t,  0)  implies  6(t)— >-0  as  t—*oo.  By  LaSalle’s  theorem,  IT  :  X  X 
Xi  x  •  •  *x  Xjv — is  a  Lyapunov  function  for  the  feedback  system,  and  the  system  is  asymptotically 
stable.  Therefore,  the  uncertain  system  is  robustly  stable.  □ 


4  Performance  Robustness  of  Uncertain  Systems 

Consider  the  following  feedback  uncertain  system  which  is  described  as  a  feedback  system  set. 


A 

G 

u 


w 


where  w  is  some  external  disturbance  vector,  and  it  is  assumed  w  £  z  is  the  regulated 

signal  vector.  The  nominal  plant  G  has  the  following  realization 

{x  =  f(x)  +  gi(x)u  +  g2(x)w 

y  =  hx(x)  +  kn(x)u  +  k12(x)w  (29) 

z  =  h2(x)  +  k2  i(x)u  +  k22(x)w 

and  the  uncertainty  structure  is  described  by  the  set 

A  :=  {A  =  block-diag{Ai,  A2,  •  •  • ,  Ajv}  :  A,-  is  nonlinear  time  varying  causal  system  for  each  *}. 

(30) 
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Note  that  the  unstructured  uncertainty  is  a  special  case  for  N  =  1.  In  the  case  of  interests  here, 
just  the  case  where  the  uncertainty  has  finite  gain  is  considered.  In  particular,  it  is  assumed  that 
all  admissible  uncertainties  are  in  the  following  set. 

BA  :=  {A  =  block-diag{Ai,  A2,  •  •  •,  Ajv}  :  A  e  A  and  has  £2-gain  <  1}  (31) 

If  A  :=  block-diag{Al5  A2,  •••,  Ajv}  €  BA,  then  each  nonlinear  system  A,-  has  £2-gain  <  1- 

Definition  4.1  The  uncertain  system  depicted  above  satisfies  robust  performance  if  for  each 
A  G  BA,  the  corresponding  feedback  system  is  well  posed  and  has  C2-gain  <  1,  i.e., 

[  (IK0ll2-IK0ll2)*<° 

Jo 

for  all  T  G  E+;  in  addition,  it  is  asymptotically  stable  around  0  for  w  —  0. 

In  this  section,  we  will  examine  under  what  conditions,  the  uncertain  system  depicted  above 
has  robust  performance. 


4.1  Case  I:  General  Nominal  Systems 


By  “general”,  we  mean  here  that  there  is  no  other  special  structural  constraint  on  nominal  systems, 
except  that  the  feedback  structures  of  the  uncertain  systems  are  required  to  be  well-posed  for  each 
uncertainty  structure  A  G  BA.  The  following  assumption  is  made. 

Assumption  4.2  For  each  A  G  BA,  A;  (i  G  {1, 2,  •  •  • ,  N})  has  the  following  realization: 

f  ii  ~  fiifi.il  £  Vi) 

(  h lifii i  £  y i ) 

which  evolves  on  X,  and  has  a  unique  asymptotically  stable  equilibrium  at  0  G  X,:  for  yt  —  0;  in 
addition,  there  is  a  C1  storage  function  Ui  such  that  Ui(^i)  <  ||ysj|2  -  1 1 ti* 1 1 2  +  fiiifii)  with  some 
negative  definite  function  :  X;-^R+. 


We  take  the  similar  scaling  treatment  for  the  nominal  system  G  to  reduce  the  conservatism 
arising  from  the  structural  constraints  of  the  uncertainty.  Define  the  scaling  matrix  set  V  as 

V  :=  {block-diag{di/,  d2I,  ■  ■  ■ ,  dNI}  :  for  each  i,  d,  G  K,  di  >  0}  (32) 

It  is  easy  to  see  that  for  each  D  G  V,  A  G  BA  if  and  only  if  D~XAD  G  BA. 

We  now  define 


g{x)-=  9\{x)  92{x)  ,/i(x):= 


hfix) 

h2(x) 


,  k{x ) 


&n(a:)  k12{x) 

k2i{x)  k22{x) 


we  have  following  theorem  about  robust  performance. 

Theorem  4.3  Under  assumption  the  uncertain  system  has  robust  performance  if  there  exist 
a  positive  definite  function  V  :  X— and  a  positive  definite  matrix  Q  G  D  such  that  the  following 
NLMI  holds: 

I j{x)f(x)  +  hT(x)Qh(x) 


M(V,Q,x):-- 


1  dV 

2  dx 


L  21 


9T(x)§ir(x)  +  kT(x)Qh(x) 


(x)g(x)  +  hT  (x)Qk(x) 
kT  (x)Qk(x)  —  Q 


<  0 


(33) 


with  Q  :~ 


Q  0 
0  I 


for  all  x  G  X  \  {0}. 
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Proof  Consider  (33),  it  is  equivalent  to  the  following  inequality, 


I  0 

0  Q-1/2 


M(V,Q,x) 


I  0 
0  Q“1/2 


| £(*)/(*)  +  hT(x)Qh(x)  \%{x)g(x)Q-^  +  hT (x)Qk(x)Q~^  ' 

^Q-1/2gT(x)^f(x)  +  Q-1/2kT(x)Qh(x)Q-1/2  Q-1^kT(x)Qk(x)Q~1/2  -  I 

Define  g(x )  g{x)Q~1/ 2,  h(a;)  :=  Q1/2h(x)  and  k(x )  :=  Qll2k{x)Q~112 .  Using  Schur  comple¬ 

ments  argument  and  the  above  inequality  is  equivalent  to  the  following  two  inequalities. 

R(x)  :=  I  —  kT (x)k(x)  >  0, 


and 


At  t 

H(V,Q,x):=  —(x)f(x)  +  hT(x)h(x)+ 

+(^^-(a;)Ka:)+^T(:c)fc(a;))(J-^T(a;)^(a:)r1(a:)(^T(a;)^-(a:)+^T(a;)Ka:))  <  °-(35) 


for  all  a:  G  X  \  {0}.  Take  V  as  defined  in  the  statement;  define  u  :=  Q1^ 
then 


n 

re 


and  y  =  Ql! 2 


y 

z 


V{x)  <  p(*)||2  -  \\m\\3  +  n{V,Q,x).  (36) 

On  the  other  hand,  notice  that  Q 1^2  G  V\  denote  Q1^2  =  block-diag{<2i I,  q2I,  ■  ■  -,qNI}.  From 
the  assumption  4.2,  for  all  A  G  BA,  there  is  a  positive  definite  function  Ui  :  X,;— >R+  for  nonlinear 
system  f&A,^-1  for  each  i  G  {1,2,  •  •  -,N}  such  that 


Mb)  <  IlftOOII2  -  ll«i(0ll2  +  Mb) 


(37) 


for  some  negative  definite  function  ipi  on  X,,  where  =  qiUi,  iji  =  qt and  £,•  is  the  state  vector  of 
A*  on  Xj.  Note  that 


Ui 

"  Vi 

U2 

h 

I 

1  £  = 

UN 

Vn 

w 

z 

Therefore, 

N  1 V 

ll^ll2  =  23  Pill2  +  IMI2  »  Pll2  =  23  Pill2  +  INI2  • 

»  =  1  i  =  1 

Next,  define  a  positive  definite  function  W  on  X  x  Xj.  X  •  •  •  x  XN  as 

N 

W(x,  6,  •  •  ■ ,  6v)  =  V{x)  +  £  Ui(Zi)  (38) 

8  =  1 


0, (34) 
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So  from  (36)  and  (37),  it  follows  that 

w(x,i i,  •••,£*)<  ll^)||2  -  \\m\f  +  n(v,Q,x)  +  Y,(.Mt)\\2  -  IIMOII3  +  Mb)) 

2  =  1 

<  IM|2  -  IM|2  +  K(V,Q,x)  +  tM6)  (39) 

2  =  1 

it  1 1 2  it  | f 2 

<  ||w||  -  ||z|| 

The  latter  inequality  implies 

jT 

/  (IMI2  _  llzll2)^  ^  9 

Jo 

for  all  T  6  R+,  i.e.,  the  feedback  system  has  £2-gain  <  1. 

Next,  we  consider  the  asymptotic  stability  for  w  =  0.  In  this  case,  (39)  becomes 

N 

W{x,t i,  M|2  +  H(V,Q,x)  +  £  &(£,-) 

2  =  1 

Thence,  if  W(x,  £i,  •  ■  ■ ,  £N)  =  0,  then  H(V,Q,x)  =  0  and  ^(&)  =  0  for  each  i;  the  former 
condition  implies  x  =  0  by  (35),  and  the  latter  implies  &  =  0  by  the  assumption.  Therefore, 
W  :  X  x  Xt  X  •  •  •  X  Xjv-*1R+  is  a  Lyapunov  function  for  the  feedback  system,  and  the  system  is 
asymptotically  stable. 

Therefore,  we  conclude  that  the  uncertain  system  satisfies  robust  performance. 

□ 


Next,  we  further  relax  the  condition  for  the  last  theorem  to  get  an  alternative  characterization 
for  the  robust  performance  of  the  given  uncertain  system. 

Assumption  4.4  Consider  the  nominal  system  G,  define  a  new  system 

(  x  =  f(x)  +  g1(x)u 
\  z  =  k2{x)  +  k2i(x)u 

The  solution  for  all  possible  u(t )  under  the  constraint  z(t)  =  0  satisfies  x(t)—> 0  as  t-+  oo. 

Remark  4.5  If  the  system  is  linear,  the  above  assumption  is  just  a  characterization  of  the  system 
with  stable  transmission  zeros. 


Theorem  4.6  Under  assumptions  4-4  an<J  4-%>  the  uncertain  system  has  robust  performance  if 
there  exist  a  positive  definite  C1  function  V  :  X— and  a  positive  definite  matrix  Q  £  V  such  that 
the  following  NLMIs  hold: 


1 j(x)f(x)  +  hT(x)Qh(x)  \%(x)g(x)  +  hT(x)Qk(x ) 

2ST(x)Sdr(x)  +  kT(x)Qh(x)  kT(x)Qk(x )  -  Q 


kT(x)Qk(x )  —  Q  <  0 


with  g(x),h(x),k(x)  defined  previously  and  Q  := 


Q  0 
0  I 


for  all  x  G  X. 


(42) 
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Proof  The  proof  is  a  combination  of  the  one  for  the  previous  theorem  and  the  one  for  theorem  2.7. 
We  just  give  an  outline  here. 

Take  V  :  X-»R+  as  given,  and  U(  :  X,:^R+  as  in  the  proof  of  the  last  theorem.  Define  a  positive 
definite  function  W  on  X  X  Xj  X  ■  •  •  X  Xjy  as 

^(*,£ir-^6v)  =  n*) +!>(&)•  (43) 

2  —  1 


It  follows  that 

w{x,  fi,  •  •  • ,  6v)  <  ikir  -  pii2 + i:  uti)  (44) 

i= 1 


<  ||U7||  —  ll^ll 

The  latter  inequality  implies 

f  (iM|2  —  llzi|2)^  >  0 

Jo 

for  all  T  G  R+,  i.e.  the  feedback  system  has  £2-gain  <  1. 

Next,  we  consider  the  asymptotic  stability  for  w  =  0.  In  this  case,  (44)  becomes 


N 


w(^,ei,---,Civ)<-Ni2+E^) 


2  =  1 


(45) 


Thence,  if  W(x,^,  •  -  • ,  £N)  —  0,  then  z  —  h2{x)x  +  k2 i(x)u  —  0  and  =  0  f°r  each  *j 

the  former  condition  implies  x(t)-+0  as  oo  by  assumption  4.4,  and  the  latter  implies  £,•  =  0  by 
assumption  4.2.  By  LaSalle’s  theorem,  W  :  X  X  Xj.  X  •  •  •  X  XN->R+  is  a  Lyapunov  function  for  the 
given  closed  loop  system,  and  the  system  is  asymptotically  stable. 

Therefore,  we  conclude  that  the  uncertain  system  has  robust  performance.  □ 


4.2  Case  II:  Strictly  Causal  Nominal  Systems 

In  this  subsection,  we  consider  the  case  where  the  nominal  system  G  is  strictly  causal  with  respect 
to  the  input  u.  In  this  case,  the  well-posedness  of  the  feedback  system  is  guaranteed  for  each 
A  G  BA.  Accordingly,  the  uncertainty  class  BA  is  enlarged. 

The  nominal  system  G  has  the  following  realization. 


X  =  f(x)  +  gi(x)u  +  g2(x)w 
y  =  h1(x)+  ki2(x)w 

z  =  h2(x)  +  k2X(x)u  +  k22(x)w 


The  uncertainty  structure  is  described  by  the  set 

BA  :=  {A  =  block-diag{Al5  A2,  •  •  •,  AN}  :  A  £  A  and  has  £2-gain  <  1} 


(46) 


(47) 
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Assumption  4.7  For  all  A  G  BA,  A;  (i  G  {1,2,  ■  •  ■,  N})  has  the  following  realization: 

{  —  h'lifl-.  C  Vi) 

which  evolves  on  X,  and  has  a  unique  asymptotically  stable  equilibrium  at  0  G  X,-  for  y,  =  0;  in 
addition,  there  is  a  C1  storage  function  U,  such  that  £/,(£,•)  <  \\yi\\  —  ||wi||  ■ 


The  scaling  matrix  set  V  is  defined  as  before 


V  {block-ding^/,  d2I,  ■  ■  -,dNI}  :  for  each  i,  di  G  R,  d,  >  0} 


It  is  easy  to  see  that  for  each  D  G  V,  A  £  BA  if  and  only  if  D  XAD  G  BA. 
We  now  define 


$(*)  :=  ShO)  9-ii.x)  J  ,h(x)  := 


hi(x) 

h2(x) 


,k(x)  := 


0  ku(x) 

k2i(x)  k-22{x)  ’ 


we  have  the  following  result  about  robust  performance  analysis. 


Theorem  4.8  Under  assumption  f.  7,  the  uncertain  system  with  nominal  plant  (46)  has  robust 
performance  if  there  exist  a  positive  definite  C1  function  V  :  X—>R  and  a  positive  definite  matrix 
Q  G  V  such  that  the  following  NLMI  holds: 

^(x)f{x)  +  hT(x)Qh(x)  §§£  {x)g{x)  +  hT{x)Qk(x)]  . 

.  +  kT(x)Qh(x)  kT(x)Qk(x)  -  Q 


with  Q 


Q  0 
0  I 


for  all  x  G  X  \  {0}. 


Proof  Consider  (48),  it  is  equivalent  to  the  following  inequalities 
R(x)  :=  I  —  kT(x)k(x)  >  0, 


and 

H(V,  Q,x)  :=  ^(x)f(x)  +  hT(x)h(x)+ 

1  f)V  „  .  1  flvT 

(x)-Qfr(x)+kT (x)h(x))  <  0.(49) 

for  all  x  G  X  \  {0}. 

Take  V  :  X-^R+  U{  :  X8--*R+  as  given.  Define  a  positive  definite  function  W  on  XxXi  X  •  •  -xX^v 
as 


w(x,zu--.,zN)  =  v(x)+jrUi(zi) 

i= 1 


(50) 
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Similar  arguments  in  the  previous  section  yield  that 


\\w\f  -\\z\f  +  n(V,Q,x)  (51) 

/  II  1 1 2  I.  m2 

<  ||'IC||  —  ||.i  || 

The  latter  inequality  implies 

/T(l|w||2  -  \\z\\2)dt  >  0  (52) 

Jo 

for  all  T  E  R+,  i.e.,  the  feedback  system  has  £2-gain  <  1. 

Next,  we  consider  the  asymptotic  stability  for  w  =  0.  In  this  case,  (51)  becomes 

Thence,  if  W(x,  £i,  •  •  • ,  (N)  —  0,  then  7i(V,Q,x)  =  0,  it  in  turn  implies  x  —  0  by  (49).  But 
x  —  0  implies  y  =  0,  therefore  £,(t)^0  as  t-> oo,  for  £;  =  /;(£;,  f,  0).  By  LaSalle’s  theorem, 

W  :  X  X  Xi  X  ■  •  •  X  Xjv~^IR+  is  a  Lyapunov  function  for  the  given  closed  loop  system,  and  the 

system  is  asymptotically  stable. 

Therefore,  we  conclude  that  the  uncertain  system  is  of  robust  performance.  □ 


Next,  we  further  relax  the  condition  for  the  last  theorem  to  get  an  alternative  characterization 
for  the  robust  performance  of  the  depicted  uncertain  system  with  nominal  system  (46). 

Assumption  4.9  Consider  the  nominal  system  G,  define  a  new  system 

f  i  =  f{x)  +  9i(x)u 

\  z  =  h2{x)  +  k2\{x)u 

The  solution  for  all  possible  u(t )  under  the  constraint  z(t)  =  0  satisfies  x{t)  =  0  for  all  t  E  R+. 

It  is  noted  that  in  the  linear  case,  the  above  assumption  corresponds  to  the  condition  that  the 
system  has  no  transmission  zero. 


Theorem  4.10  Under  assumptions  f.9  and  f.7,  the  uncertain  system  has  robust  performance  if 
there  exist  a  positive  definite  C1  function  V  :  X-h>R  and  a  positive  definite  matrix  Q  E  T>  such  that 
the  following  NLMIs  hold: 

frO  )/(*)  +  hT(x)Qh(x)  ^(x)g(x)  +  hT(x)Qk(x) 

.  1 9T{x)^jr(x)  +  k'r (x)QKx)  kT(x)Qk(x)  -  Q 

kT(x)Qk(x )  -  Q  <  0 


(53) 

(54) 


with  B(x),  C(x),  D(x)  defined  previously  and  Q  := 


Q  0 
0  I 


for  all  i£X, 
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Proof  The  proof  is  a  combination  of  the  one  for  previous  theorem  and  the  one  for  theorem  2.7. 
We  just  give  an  outline  here. 

Take  V  :  X-»R+  as  given,  and  £/,-  :  X2  ^E+  as  in  the  proof  of  the  last  theorem.  Define  a  positive 
definite  function  W  on  X  X  Xt  X  ■  ■  ■  X  X,v  as 

W(z,{u  •••,6v)  =  V(x)  +  J2  U^)  (55) 

2  —  1 

It  follows  that 

w(x,6,-^6v)<IMI2-IMI2 

The  latter  inequality  implies 

rp 

[  (IM|2  —  lkl|2)^  >  0  (56) 

Jo 

for  all  T  £  R+,  i.e.,  the  feedback  system  has  £2- gain  <  1. 

Next,  we  consider  the  asymptotic  stability  for  w  =  0.  In  this  case, 

W(z,6, ••■,£*)  <  -  INI2 

Thence,  if  W(x, 6,  ■  ■  ■ , 6v)  =  0,  then  2  =  h2{x)  +  k2X(x)u  =  0  which  implies  x(t)  =  0  by 
assumption  4.9.  But  x  =  0  implies  y  =  0,  therefore  &(i)— »•  0  as  t-* oo,  for  £  =  t,  0).  By 
LaSalle’s  theorem,  W  :  X  x  Xt  x  •  •  •  X  XN^R+  is  a  Lyapunov  function  for  the  given  closed  loop 
system,  and  the  system  is  asymptotically  stable. 

Therefore,  the  uncertain  system  satisfies  robust  performance.  □ 


5  On  Robustness  Synthesis 


In  the  robustness  analysis  results  in  the  last  two  sections,  the  robustness  conditions  of  uncertain 
systems  are  essentially  the  small  gain  conditions  for  (scaled)  nominal  systems  modulo  some  ap¬ 
propriate  stabilizing  conditions.  So  the  robustness  synthesis  can  be  pursued  by  combining  the 
robustness  analysis  results  in  the  last  two  sections  with  the  treatments  of  Hoo -control  synthesis  (see 
for  example  [31,  3,  18]).  We  just  take  the  performances  robustness  synthesis  problem  as  an  example, 
the  other  problems  can  be  done  similarly.  Technically,  we  closely  follow  the  treatments  in  [18,  17], 
so  we  just  take  the  state  feedback  case  as  an  example.  The  output  feedback  case  can  be  done 
similarly  by  just  modifying  the  treatments  in  [18,  17].  It  is  noticed  that  the  robust  stabilization 
with  unstructured  uncertainty  is  also  considered  in  [32]. 

Consider  the  following  feedback  uncertain  system  which  is  described  as  a  feedback  system  set. 


u 


w 

uc 
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where  w  G  £2(R+)  is  some  external  disturbance  vector,  2  is  the  regulated  signal  vector,  ym  is  the 
measured  output  vector,  and  based  on  which  the  control  input  vector  ue  is  produced.  The  nominal 
plant  G  has  the  following  realization 

x  =  f(x)  +  gi(x)u  +  g2{x)w  +  g3(x)uc 

Q  .  V  =  h^x)  +  ku(x)u  +  k12(x)w  +  k13(x)uc  , 

'  j  2  =  h2(x)  +  k21(x)u  +  k22(x)w  +  k23(x)uc 
Vm  =  h3(x )  +  k31(x)u  +  k32(x)  w  +  k33(x)uc 

where  G  C°,  and  /( 0)  =  0,h,-(0)  =  0,  for  i,j  =  1,2,3.  In  this  section,  the  state  vector 

of  the  nominal  system  is  directly  measured,  i.e.,  ym  =  x;  the  uncertainty  structure  is  described  by 
the  set 

A  :=  {A  =  block-diag{Ai,  A2,  •  •  •,  Ajv}  :  A,-  is  nonlinear  time  varying  causal  system  for  each  i). 
And  the  admissible  uncertainties  are  in  the  following  set. 

BA  :=  {A  =  block-diag{Ai,  A2,  •  •  •,  Ajv}  :  A  G  A  and  has  £2-gain  <  1}  (58) 

The  performance  robustness  synthesis  problem  by  state  feedback  is  defined  as  follows. 

Definition  5.1  (State  Feedback  Synthesis  Problem)  Find  a  state  feedback  law  uc  =  K(x) 
with  K  G  C°  and  K( 0)  =  0  for  the  uncertain  system  depicted  above  such  that  the  closed  loop 
uncertain  system  satisfies  robust  performance. 

If  uc  =  K( x)  is  a  state  feedback  law,  then  the  closed  loop  uncertain  system  is  as  follows. 


with 

x  =  (/(a?)  +  g3(x)K(x))  +  gi(x)u  +  g2(x)w 
Gf  :  y  =  (hi(x)  +  k13(x)K(x))  +  kn(x)u  +  ki2(x)w 
z  =  (h2(x)  +  k23(x)K(x))  +  k2i{x)u  +  k22(x)w 

Define  the  scaling  matrix  set  V  as 

V  :=  {block-diag{di/,  d2I ,  •  •  • ,  dNI}  :  for  each  i,  dt  G  R,  dt  >  0}  (59) 

It  is  easy  to  see  that  for  each  D  G  V,  DA  =  AD  for  A  G  A;  and  A  G  BA  if  and  only  if 
D_1AD  G  BA. 

Next,  we  consider  two  cases  about  robustness  synthesis  by  state  feedback. 
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5.1  State  Feedaback  Solutions 

Consider  the  uncertain  system  with  the  nominal  plant  as  (57).  Define 


g(x)  :=  §i(x)  g2(x)  ,h(x):= 


hi(x) 
h2(x)  ’ 


kx(x)  := 


kn(x)  k12(x) 
k21(x)  k22(x) 


,k2(x) 


kia(x) 

k23(x) 


The  following  structural  constraints  are  imposed. 


Assumption  5.2  kfix)  =  0,  and  klf(x) 
x  e  x. 


h(x)  k2(x) 


0  Ro(x)  where  R0(x )  >  0  for  all 


It  is  noted  that  if  A  :=  block-diag{ A1?  A2,  •  •  • ,  Ajv}  E  BA,  then  each  nonlinear  system  A*  has 
£2-gain  <  1.  Furthermore,  we  have  the  following  stronger  assumption. 

Assumption  5.3  For  each  A  6  BA,  A,-  (i  6  {1,2,  •  ■  -,N})  has  the  following  realization: 

f  ti  = 

\  Ui  =  hi(fi,  t,  yf) 

which  evolves  on  X<  and  has  a  unique  asymptotically  stable  equilibrium  at  0  G  X,:  for  —  0;  in 
addition,  there  is  a  C1  storage  function  U ,•  such  that  Ufifi)  <  \\yi\f  -  U-m* ||2 . 

We  first  have  the  following  lemma. 

Lemma  5.4  Consider  the  system  defined  in  (51)  with  the  structural  assumption  5.2.  The  following 
two  statement  are  equivalent. 

(i)  There  exist  a  C°  vector-valued  function  K(x)  on  X,  a  C1  positive  definite  function  V  : 
X— >R+,  and  a  positive  definite  matrix  Q  &V  such  that  the  following  NLMI  holds, 

j^(x)(f(x)  +  g3(x)K(x))  +  (h(x)  +  k2(x)I((x))TQ(h(x)  +  k2(x)K(x))  j(x)g(x) 

[  bT  -Q  J 

for  all  x  G  X  \  {0}. 

(ii)  There  exist  a  C1  positive  definite  function  V  :  X— >R+  and  a  positive  definite  matrix  Q  £  V 
such  that  the  following  NLMI  holds, 

dV  1  dV  F)VT 

-Q-(x)f(x)+  --^(x)(g(x)Q  1gI  (x)  -  g3(x)R-1(x)gf(x))-^-- (x)  +  hT (x)Qh(x)  <  0. 
for  all  x  E  X  \  {0}. 

Moreover,  if  (ii)  is  true,  then  a  state  feedback  function  K(x)  makes  (i)  true  is  as  follows. 

K(x)  =  -l~R-\x)gl(x)^~{x). 
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Proof  Note  that  the  NLMI  in  statement  (i)  is  equivalent  to  the  following  Hamilton-Jacobi  in 
equality, 

^rOX/OO  +  g3(x)K(x))  +  x)g(x)Q~1gT(x)^-(x)+ 

+(h(x)  +  k2(x)K(x))TQ(h{x)  +  k2(x)K(x))  <  0 

for  all  x  E  X  \  {0}.  By  the  same  arguments  as  in  [17,  Theorem  4.1],  the  conclusion  follows.  □ 
The  main  result  in  this  subsection  is  stated  as  follows. 


Theorem  5.5  Consider  the  uncertain  system  with  nominal  plant  as  (57).  Under  assumptions  5.3 
and  5.2,  the  state  feedback  robust  performance  synthesis  problem  has  a  solution  if  there  exist  a 
positive  definite  C1  positive  definite  function  V  :  X— >M+  and  a  positive  definite  matrix  Q  £  V  such 
that  the  following  NLMI  holds, 

^(x)f(x)  +  ~^(x)(g(x)Q~1gT(x)  -  g3{x)R3l(x)gl{x))^-{x)  +  hT (x)Qh(x)  <  0.  (60) 

for  all  x  €  X\{0}.  Moreover,  if  (V(x),Q)  is  such  a  pair  of  solutions,  then  a  state  feedback  function 
K(x)  makes  the  closed  loop  system  has  a  robust  performance  is  K(x )  =  —T,Rf1{x)gf{x)^(x). 


Proof  Let  (V(x),Q)  be  as  in  the  theorem.  By  the  preceding  lemma,  there  exist  a  C°  matrix 
valued  function  K(x)  on  X  defined  as  K(x)  —  —  ^RQ1(x)gJ(x)^-(x),  such  that 


%{x){f{x)  +  9z{x)K{x))  +  (h(x)  +  k2(x)K(x))TQ(h(x)  +  k2(x)K(x)) 

|  gT(x)B£(x) 


1  dV 

2  dx 


(x)g(x) 

-Q 


<  0(61) 


for  all  x  E  X  \  {0}.  On  the  other  hand,  take  uc  =  K(x)  as  a  state  feedback  law,  so  the  closed  loop 
nominal  system  is  as  follows. 


Gp  ■  < 


x  =  (f(x)  +  g3(x)K(x))  +  gi(x)u  +  g2(x)w 
V  =  (hfix)  +  k13(x)K(x)) 
z  =  (h2{x)  +  k23(x)K(x)) 


By  theorem  4.8,  the  closed  loop  uncertain  system  satisfies  robust  performance. 


□ 


Note  that  the  above  characterization  is  not  convex  in  general.  In  the  next  subsection,  we  will 
give  a  convex  characterization  which  have  some  computationally  appealing  property. 


5.2  A  Convex  Characterization  for  State  Feedback  Solutions 


In  this  section,  in  stead  of  the  nominal  plant  (57),  the  following  nominal  plant  is  examined. 


x  =  A(x)x  +  B\{x)u  +  B2{x)w  +  B3(x)uc 
y  =  Ci(x)x  +  Dn(x)u  +  Dl2{x)w  +  Dl3(x)uc 
z  =  C2{x)x  +  D21(x)u  +  D22(x)w  +  D23(x)uc 
.  Vm  =  X 


(62) 


where  A,  B,,Cj,  D,j  are  C°  matrix- valued  functions. 
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We  now  define 


B(x)  := 

Bfix) 

B2(x) 

,C(x)  := 

'  Cfix)  ' 
Cfix) 

,D(x)  := 

Dn(x)  Dl2(x) 
D2i(x)  D22(x) 

B(x) := 

Bl(x) 

DU*) 

DU*) 

and  let  A f(B(x))  be  the  distribution  on  X  which  annihilates  all  of  the  row  vectors  of  B(x). 

It  is  noted  that  if  A  :=  block-diag{Ai,  A2,  •  •  • ,  XN}  G  BA,  then  each  nonlinear  system  A,-  has 
£2-gain  <  1.  Furthermore,  we  have  the  following  stronger  assumption  on  the  uncertainty. 

Assumption  5.6  For  each  A  G  BA,  A;  (i  G  {1,2,  •••,  N})  has  the  following  realization: 

f  F  —  fi  ( s  /  *  l  *  Hi ) 

(  —  hi(£i,t,  yi) 

which  evolves  on  X,  and  has  a  unique  asymptotically  stable  equilibrium  at  0  G  X*  for  yt  =  0;  in 
addition,  there  is  a  C1  storage  function  Ui  such  that  Ufifi)  <  ||?/i||2  —  1 1 1 1 2  +  i’i(fii)  with  some 
negative  definite  function  tfi  :  X;— >R+. 


We  first  have  the  following  lemma. 


Lemma  5.T  The  following  two  statements  are  equivalent. 

(i)  There  exist  a  C°  matrix  valued  function  F(x),  a  positive  definite  matrix-valued  function 
P  :  X-^R”*",  and  a  positive  definite  matrix  Q  G  V  such  that  the  following  NLMI  holds: 


Atf{x)P{x)  +  P{x)Af{x)  +  Cf(x)QCf 
BT(x)P(x)  +  Dt{x)QCf{x) 


(*)  P(x)B(x)  +  Cf{x)QD(x ) 
DT{x)QD{x )  -  Q 


<  0 


with  Q 


Q  0 
0  I 


for  all  x  G  X,  where 


AF(x)  =  A(x)  +  B3(x)F(x ),  CF(x) 


C^x)  +  D13(x)F(x) 
C2(x)  +  D23(x)F{x) 


(63) 


(ii)  There  exist  a  positive  definite  matrix- valued  function  X 
matrix  Y  G  D  such  that  the  following  NLMI  holds: 


X- 


<n  and  a  positive  definite 


B T 


0) 


A(x)X(x)  +  X{x)AT{x)  +  BT(x)YB(x) 
C(x)X(x)  +  D(x)YBT(x) 


X(x)CT(x)  +  B{x)YDt{x) 
D(x)YDt(x)-Y 


B±(x)  <  0 
(64) 


with  Y  := 

A f(B(x))  for  all  x  G  X. 

Moreover,  If  any  one  of  the  above  statments  holds,  then  the  solutions  of  the  other  NLMI  can  be 
chosen  such  that  P(x)  =  X~1(x)  and  Q  =  Y~l. 


and  B±(x)  is  a  C°  matrix-valued  function  on  X  such  that  span(B±(x ))  = 
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Proof  Use  the  similar  argument  in  [18,  Section  4]. 


□ 


It  is  noted  that  the  NLMI  (64)  is  affine  in  unknown  P(x)  and  Q.  We  have  following  theorem 
which  convexly  characterizes  robust  performance  synthesis  by  state  feedback. 

Theorem  5.8  Consider  the  uncertain  system  with  nominal  plant  defined  as  (62).  Under  assump¬ 
tion  5.6,  the  state  feedback  robust  performance  synthesis  problem  has  a  solution  if  there  exist  a 
positive  definite  matrix-valued  function  X  :  X-+Rnxn  and  a  positive  definite  matrix  Y  £  V  such 
that  the  NLMI  (64)  holds  for  all  x  £  X,  and  §7(2)  =  2xTX~1(x)  for  some  C1  function  V  on  X 
with  P(0)  =  0. 


Proof  Let  (X(x),Y)  be  as  in  the  theorem.  By  the  preceding  lemma,  there  exist  a  C°  matrix 
valued  function  F{x)  on  X,  such  that 


A^(x)X  1(x)  +  X  l(x)AF(x)  +  CF 

BT(x)X-\x)  + Dt(x)Y~1Cf(x) 


f(x)Y 


-1 


CF(x)  X-\x)B{x)YCtf(x)Y-1D(x) 
DT(x)Y-1D(x)-Y~1 


<  0 


(65) 


On  the  other  hand,  take  uc  =  F(x)x  as  a  state  feedback  law,  so  the  closed  loop  nominal  system 
is  as  follows. 


Gf  '■  < 


x  =  {A{x)  +  B3(x)F(x))x  +  Bi(x)u  +  B2(x)w 
y  =  (Ci(a:)  +  Di3(x)F{x))x  +  Dn(x)u  +  D12(x)w 
z  =  ( C2(x )  +  D23(x)F(x))x  +  D21(x)u  +  D22(x)w 


We  now  claim  the  closed  loop  uncertain  system  satisfies  robust  performance.  In  fact,  take 
(. P(x),Q )  =  (W_1(a;),y-1).  By  the  remarks  in  section  2.3,  (65)  implies  that 


xT(A£(x)P(x)  +  P(x)Af(x)  +  CJ.(x)QCf(x))x  xt{P{x)B(x)  +  CTF{x)QD{x )) 
( BT(x)P(x )  +  Dt(x)QCf{x))x  Dt(x)QD(x)  -  Q 


<  0 


The  conclusion  therefore  follows  from  theorem  4.3. 


□ 


6  Computational  Issues 


We  address  computational  issues  for  robustness  analysis  and  synthesis  in  this  section.  From  the 
development  of  the  theory,  it  is  noted  that  the  computation  about  robustness  analysis  and  synthesis 
involves  solving  some  NLMIs.  To  be  more  concrete,  we  take  the  following  NLMI  with  respect  to 
(V(x),Q)  (i.e.  (21)  or  (33))  as  an  example, 


M(V,Q,x):= 


f j(x)f(x)  +  hT(x)Qh(x)  lSV(x)g(x)  +  hT(x)Qk(x) 

29T(x)S£-(x)  +  kT(x)Qh(x )  kT(x)Qk(x)  -  Q 


for  all  2:  6  X.  The  computation  procedure  is  therefore  divided  into  the  following  two  steps, 


e  Find  ( p(x),Q )  for  some  C°  vector  valued  function  p  :  X— dR”  with  p(0)  =  0  and  positive 
definite  matrix  such  that 


2  PT(x)f(x)  +  hT(x)Qh(x)  pT(x)g(x)  +  hT(x)Qk(x) 
gT(x)p(x)  +  kT  (x)Qh(x)  kT(x)Qk(x )  —  Q 
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•  Check  if  there  is  a  Lyapunov  function  V  :  X^R+  such  that  =  2 pT(x)  for  iEX. 

Next,  we  pursue  these  two  issues.  The  reader  is  also  referred  to  [18]  for  a  detailed  discussion 
about  a  special  case.  The  technique  used  here  is  the  same  as  that  in  [18]. 


6.1  Solutions  of  NLMIs 


We  first  consider  the  solution  about  (66).  As  discussed  in  section  2.3,  if  we  assume  f(x)  = 
A(x)x,g(x)  =  B(x),h(x)  =  C(x)x,  k(x)  =  D(x),  and  p(x)  =  P(x)x  in  (66),  then  (66)  is  im¬ 
plied  by  the  following  NLMI  (which  is  a  more  conservative  characterization  about  robustness), 


M(P,Q,x) 


AT(x)P(x)  +  PT(x)A(x)  +  CT(x)QC(x)  PT(x)B{x )  +  CT{x)QD(x) 
BT{x)P{x)  +  DT(x)QC(x)  Dt(x)QD(x)  -  Q 


<  0. 


(67) 


In  this  subsection,  we  need  to  find  a  C°  matrix-valued  function  P  :  X— >RnX"  and  a  positive 
definite  matrix  Q  such  that  M(P,Q,x)  <  0.  In  other  word,  if  (P(x),  Q )  is  a  solution  to  (67),  then 
( P(x)x ,  Q )  is  a  solution  to  (66). 

Under  some  regularity  conditions,  the  existence  of  continuous  and  positive  definite  solutions  to 
NLMIs  is  also  justified.  Given  a  matrix  valued  function  S  :  X— »Rnxn  with  S(x)  >  0  for  all  x  £  X. 
Let  R(x)  Q  —  DT(x)QD(x )  >  0,  define  a  state-dependent  Hamiltonian  H  :  X^R2nx2n  as 


H(x)  := 


A(x)  0 

CT(x)QC(x)  -  S(x)  - AT(x ) 


+ 


+ 


B(x) 

- CT(x)QD(x ) 


R  x(a;)  DT(x)QC(x)  BT(x )  . 


The  following  result  is  essentially  from  [13,  lemma  2.4]. 


(68) 


Theorem  6.1  A4(P,Q,x)  <  0  has  non-negative  definite  solutions  P(x)  >  0  and  Q  >  0  if  and  only 
if  the  state- dependent  Hamiltonian  H  :  X-*R2nx2n  defined  in  (68)  for  some  matrix-valued  function 
S  :  X^R"xn  with  S(x)  >  0  for  all  x  G  X  is  in  dom(Ric),  i.e.  H(x)  6  dom(Ric)  for  each  x  6  X. 
Moreover,  P(x )  :=  Ric(H(x ))  >  0  is  such  a  solution  with  Q  >  0  as  given.  In  addition,  if  for  each 

x  e  x, 


n  — 1 


P|  ker {C(x)Ai{x))  =  0,  C{x) 


i= 0 


Q1/2C(x) 

S^(x) 


this  solution  is  positive  definite,  i.e.,  Ric(H(x))  >  0. 


The  above  theorem  implies  that  under  the  condition  H(x)  6  dom(Ric)  for  each  x  G  X,  the 
NLMI  M(P,Q,x)  <  0  has  non- negative  definite  solutions  P(x )  >  0  and  Q  >  0.  The  following 
theorem  further  shows  that  such  solutions  can  be  chosen  to  be  continuous  in  the  case  of  interest  in 
this  paper. 


Theorem  6.2  Suppose  the  matrix  inequality  M(P,Q,x)  <  0  has  a  positive  definite  solution  P(x) 
for  each  x  G  X  and  Q  >  0,  then  there  exists  a  C°  matrix-valued  function  P  :  X— >R”X"  with 
P(x)  —  PT(x)  >  0,  such  that  M(P(x),Q,x)  <  0  for  all  x  G  X. 
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Proof  Consider  the  NLMI: 

M(P,Q,x)<  0. 

with  x  E  X,  where  M  :  RnXn  x  Rpxp  x  X^S(R(n+p:)x(n+p:))  is  continuous  and  satisfies 

N  N  N 

M(J2akpk,^2oikQk,x)  =  J2akM(Pk,Qk,x)  (69) 

k= 1  *=1  k= 1 

for  all  ak  >  0  with  i  a> *  =  1- 

By  assumption,  there  exists  a  positive  definite  matrix  Q  E  Rpxp,  for  each  x  E  X,  there  is  a 
positive  definite  Px  E  RnXn  such  that 

M(Px,Q,x)  <  0. 

By  continuity  of  M  with  respect  to  x,  there  is  a  rx  >  0  such  that  for  all  x0  E  N(x)  :=  {x0  : 
||*o  -  *||  <  rx}, 

M(Px,Q,x  o)<0.  (70) 

On  the  other  hand,  {N(a;)}|x€X  is  an  open  covering  of  X ,  i.e., 

^  C  1J  N(x)  (71) 

*6X 

Since  the  space  1"  is  paracompact,  there  is  a  locally  finite  open  subcovering  {N,}  |ieI  for  some  index 
set  I  which  refines  {N(x)}|x€X.  By  (70),  Pi  G  Rnxn  is  taken  to  be  positive  definite  for  each  i  G  I 
such  that 

M(Pi,Q,x)<  0.  (72) 

for  all  x  G  N,-. 

It  is  known  by  the  standard  results  of  continuous  partitions  of  unity  that  there  is  a  locally 
Lipschitzean  partition  of  unity  {V’i}|i€i  to  X  subordinated  to  the  covering  {N*-} |iGi ;  i.e.,  fa  is  locally 
Lipschitzean  and  non-negative  with  support  Supp(ipi)  C  N;  for  each  i  G  I,  and 

£>*(*)  =  l,Vz  E  X.  (73) 

i€  I 

Define  a  matrix-valued  function  P  :  X-^P(Rnxn)  as 

P(*)  =  X>(*)fl,V*GX,  (74) 

»e  i 

which  is  positive  definite  and  continuous  since  it  is  locally  a  finite  sum  of  continuous  positive  definite 
matrix- valued  functions. 

It  follows  from  (73),  (74)  and  (69)  that 

M(P(x),Q,x)  =  M(£fa(x)Pt,Q,x)  =  Yl^x^M(Pi,Q,x)  <  0 
igl  iel 

The  last  equality  holds  since  the  sum  is  finite  for  each  x  G  X. 

Thence,  the  constructed  C°  matrix-valued  function  P  :  X-+P(R"xn)  in  (74)  is  positive  definite 
and  is  a  solution  to  A 4(P(x),Q,x)  <0.  □ 
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Remark  6.3  The  similar  technique  can  be  used  to  examine  the  continuous  solutions  to  the  follow¬ 
ing  NLMI: 

xT(AT(x)P(x)  +  PT(x)A(x)  +  CT(x)QC(x))x  xt{Pt(x)B(x)  +  CT(x)QD(x))  ,  , 

(BT(x)P{x)  +  DT(x)QC{x))x  DT(x)QD(x)-Q  J  ' 

for  all  x  G  X  \  {0}.  As  a  matter  of  fact,  we  can  get  a  continuous  solution  P(x)  onX.\  {0}  to  the 
above  NLMI  using  the  similar  arguments  as  in  the  preceding  proof.  If  limE_*0  P{x)  =  P0  £  Rnxn, 
then  the  matrix  valued  function  P(x)  on  X  \  {0}  can  be  continuously  extended  to  X  by  defining 
P( 0)  =  P0.  The  extension  is  a  solution  to  (75)  on  X. 

A  nice  convex  property  for  NLMIs  is  stated  by  the  following  proposition  whose  proof  is  easy 
and  omitted  here. 

Theorem  6.4  The  C°  solutions  ( P(x),Q )  to  NLMI  M(P,Q,x)  <  0  such  that  P  :  X— >R"xn  and 
Q  >  0  form  a  convex  set;  the  subset  of  solutions  (P(x),  Q)  such  that  P(x)  is  C°  non-negative  definite 
with  §7(2)  =  2 xT P(x)  for  some  function  V  :  X^R  is  convex;  the  subset  of  solutions  ( P(x),Q ) 
such  that  P(x)  is  C°  positive  definite  with  §7(2:)  =  2 xT P(x)  for  some  function  V  :  X-+R  is  also 
convex. 


6.2  Existence  of  Lyapunov  Functions 

As  mentioned  earlier,  the  existence  of  solution  to  NLMI  (66)  is  not  enough  to  give  positive  assertion 
about  £2-performance,  robust  stability,  and  robust  performance;  some  additional  requirement  is 
imposed  in  this  paper,  i.e.  there  is  a  C1  Lyapunov  function  V  :  X— dR.+,  such  that 

9V<  ^  o  77  X 

for  all  x  G  X.  In  this  subsection,  we  will  examine  explicitly  when  it  is  the  case  for  a  class  of  such 
solutions. 

The  following  result  is  quite  standard,  the  reader  is  referred  to  [5]  and  [18]  for  the  proofs. 


Proposition  6.5  Suppose  a  vector-valued  function  p  :  X— >Rn  is  of  class  C1 ;  letp(x)  —  (pi(x),  ■  ■  -,pn(x)]T 
for  i£X.  Then  there  exists  V  :  X— >R  such  that 

dV,  .  „  T/  , 

&  = 2p  ^ 
if  and  only  if 


(76) 


for  all  x  G  X  and  i,j  =  1,2,  •••,«,  Moreover,  if  (76)  holds,  then  an  function  V  :  X-+R  with 
E(0)  =  0  is  given  by 


V(x)  =  2xt 


(77) 


In  addition,  if  p(x)  =  P(x)x  for  some  positive  definite  matrix-valued  function  such  that  P(x),  then 
V(x)  is  also  positive  definite  function. 
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For  a  class  of  solutions  constructed  in  the  proof  in  theorem  6.2,  we  can  similarly  characterize  the 
existence  of  the  Lyapunov  function.  The  matrix-valued  function  P  :  X->P(R"X"),  which  satisfies 
JH(P,x)  <  0,  is  constructed  as  (74) 

P(x )  =  J2Mx)Pi,Vx  E  X, 

*ei 

for  some  index  set  I,  where  is  a  partition  of  unity  of  X  and  Pi  =  Ff  >  0.  Notice  that  the 

above  summation  is  locally  finite.  Similar  argument  leads  to  the  following  theorem. 


Theorem  6.6  Suppose  the  matrix  valued  function  P  :  X^P(R"X")  defined  by 


P(x )  =  J2Mx)Pi 

ie  i 


(78) 


with  % 4  :  X^R+  being  of  class  C1  and  P{  E  P(R”xn)  for  i  El  satisfies:  A4(P,  x)  <  0  for  all  x  E  X; 
let  Vi(x')  =  xT PiX  for  all  i  E  I.  There  exists  a  C2  function  V  :  X^R  such  that  |j(»)  =  2 xTP(x) 
if  and  only  if 


E 


dxt 


w  =  E 

*ei 


dfii 

dx, 


(x) 


(79) 


for  all  x  E  X  and  j,l  E  {1,2  ,■■■  ,n}  with  j  /  /. 


Notice  that  the  summation  in  (79)  is  finite  for  each  *  E  X. 


Proof  Consult  [18]. 


□ 


6.3  Further  Remarks  and  An  Example 

The  above  treatments  about  robustness  analysis  and  synthesis  are  in  terms  of  NLMIs,  which  are 
pointwise  LMIs  on  state  set  X,  modulo  some  additional  constraints  on  the  solutions.  From  the 
proof  of  Theorem  6.2,  we  know  that  if  X  is  bounded,  then  we  only  need  to  solve  a  finite  number  of 
LMIs  to  get  the  solution  for  the  NLMI  on  X.  This  is  different  from  some  other  recent  treatments 
as  in  [6]. 

In  the  light  of  the  notion  of  global  linearization  of  nonlinear  systems  developed  by  Liu  et  al  [15], 
in  [6],  the  authors  view  the  considered  nonlinear  system  as  a  parametric  uncertain  system,  i.e.,  in 
the  current  case,  the  coefficient  matrices  in  (67)  are  assumed  in  a  convex  set: 


[A(x),B(x),C(x),D(x)}£Co{[Ai,Bi,Ci,Di]\iHli2i..,L}},Vx£X, 

where  Co  stands  for  the  convex  hull.  In  this  case,  a  constant  solution  ( P,Q )  E  Rnx"  X  R4X?  to 
(67)  is  sought.  Therefore,  if  for  allf  E  {1,  2,  •  •  • ,  L} 

'  AjP  +  Pi:  A i  +  CjQCi  PTBi  +  CjQDi  ] 

BJP  +  DjQCi  DjQDi  -Q  ~  U‘ 


have  a  common  constant  solution  ( P,Q ).  Then  ( P,Q )  is  also  a  solution  to  (67),  i.e., 


AT(x)P  +  PTA(x)  +  CT(x)QC(x)  PtB(x )  +  CT(x)QD(x ) 
BT{x)P  +  DT(x)QC(x)  DT(x)QD(x)-Q 


<  0. 
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The  solution  automatically  satisfies  the  condition  (76),  and  the  corresponding  Lyapunov  function 
is  V (x)  =  xT  Px. 

This  treatment  suggests  a  tractable  algorithm  to  get  local  solutions,  which  can  be  used  to  seek 
constant  solutions  on  each  partitioned  state  set  N,-  in  the  proof  of  theorem  6.2.  However,  this 
approach  generally  leads  to  conservative  results  if  the  prescribed  state  set  is  large  enough.  This 
can  be  seen  in  the  following  example. 

Example  6.7  Consider  an  uncertain  feedback  system  with  diagram  as  follows. 


Where  P  is  the  nonlinear  plant;  K  is  the  controller  such  that  the  output  z  is  supposed  to  be 
regulated;  ym  is  the  measured  output,  based  on  wdiich  the  control  action  uc  is  produced;  w  is  the 
disturbance  from  the  actuator;  and  u  is  the  disturbance  from  the  sensor  which  is  generated  by 
u  =  A uc  with  A  being  the  bounded  causal  uncertainties.  The  robustness  analysis  problem  is 
to  check  that  for  a  given  controller  K,  whether  or  not  the  influence  of  the  noises  w  on  the  regulated 
output  z  is  reduced  to  the  required  degree  for  all  possible  A. 

To  formulate  this  problem,  all  the  signals  are  considered  in  space  £2[0,  oo).  We  will  check  that 
given  7  >  0  and  K,  does 

f  \\z\\2  dt  <  72  [  ||w|j2  dt,  VT  6  M+, 

Jo  Jo 

for  all  admissible  A  ? 

In  this  example  ,  the  plant  has  the  following  realization: 

{x  =  ex(u  +  ue) 

Z  =  X  +  w 
ym  =  x  +  w 

and  the  controller  K  =  —  1,  each  admissible  uncertainty  A  has  £2-gain  <  We  will  check  if  the 
£2~gain  <  To  this  end,  the  standardized  block  diagram  for  the  closed  system  is  redrawn  as 
follows 
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with 


G  :  < 


x  =  —  exx  —  exu  +  exw 

y  =  -j=2x~7 2U 


[  Zr  =  72X+^U 

Therefore,  zr  =  \/2z.  In  this  case,  the  admissible  uncertainty  As  =  y/2A  is  assumed  to  have  the 
following  realization: 


(  £  = 

\  u  =  h(£,t,y) 


which  evolves  on  for  some  integer  /,  and  £  =  /(£,£  0)  is  asymptotically  stable  at  0  G  K;;  in 
addition,  there  exists  a  C1  storage  function  U  such  that  U(£)  <  ||?/||2  —  ||u||2  +  ip(£)  with  some 
negative  definite  function  ip  :  R(— >R+. 

It  is  sufficient  to  check  if  the  above  feedback  system  has  robust  performance.  We  first  consider 
zN  :=  ezr  =  c(^a:  +  ^=m)  for  some  positive  e  <  1.  Consider  the  following  NLMI, 


-2exP(x)  +  |(Q  +  e2)  -exP(x)  +  |(<3  +  e2)  exP(x) 
—exP(x)  +  |(Q  +  e2)  l(e2~Q)  0 

exP(x)  0  -1 


<  0 


which  implies  the  NLMI  in  theorem  4.3  and  is  equivalent  to 


-2  exP(x)  +  |(Q  +  c2)  +  e2xP\x)  -ex  P{x)  +  +  e2) 

—exP(x)  +  |(Q  +  e2)  He2~Q) 


<  0 


or 


e2  -  Q  <  0  (80) 

-(2  +  Q-  e2)e2xP2(x)  +  4Q2exP(x)  -  2 Q(e2  +  1)  >  0.  (81) 

There  exist  positive  solutions  ( Q,P(x ))  to  the  above  two  inequalities;  they  satisfy  Q  —  1  and 
j^e~x  <  P{x)  <  e~x.  Hence, 

f  \\zn\\2  dt  <  (  ||ru||2  dt,VT  G  R+, 

Jo  Jo 

Therefore,  the  £2-gain  for  the  closed  loop  system  <  ^  for  all  e  <  1,  which  in  turn  implies  the 
£2-gain  <  ^=.  However,  it  is  obvious  that  there  is  no  constant  pair  ( P,Q )  which  satisfies  both  (80) 
and  (81)  for  all  x  G  R. 

7  Concluding  Remarks:  Robust  BIBO  Stability  and  the  Scaled 
Small  Gain  Theorem 

This  paper  deals  with  the  robustness  analysis  and  synthesis  for  nonlinear  systems  in  the  state 
space.  The  characterizations  are  in  terms  of  NLMIs  which  offer  some  computationally  attractive 
properties.  In  an  input-output  point  of  view,  the  characterizations  imply  that  both  (structured) 
uncertainty  and  a  (scaled)  nominal  plant  have  £2-gain  <1.  It  is  noted  that  the  scaling  treatment 
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for  the  robust  performance  problem  or  robust  stability  analysis  in  the  structured  uncertainty  case 
is  a  reasonable  way  to  reduce  the  possible  conservatism  arising  from  the  structural  constraints  of 
the  uncertainty. 

As  a  matter  of  fact,  the  robust  stability  and  performance  conditions  in  this  paper  can  be 
strengthened  using  BIBO  stability  instead  of  asymptotic  stability  and  a  scaled  small  gain  theorem. 
Under  weaker  assumptions  on  the  uncertainty  the  scaled  small  gain  theorem  is  necessary  as  well  as 
sufficient,  and  can  be  related  directly  to  conditions  of  Theorems  3.5  and  4.3.  To  examine  the  robust 
stability  of  a  uncertain  system  in  the  input-output  setting,  consider  the  following  block  diagram, 


u  2 


where  the  uncertainty  A  is  causal  and  possibly  nonlinear  and  time- varying  with  £-gain  <  1,  the 
nominal  system  G  is  nonlinear  time-invariant  and  causal.  For  the  purposes  of  this  section,  no 
realizations  for  A  or  G  need  be  assumed.  The  feedback  structure  is  assumed  to  be  well-posed  for 
any  admissible  uncertainty.  The  robust  BIBO  stability  requires  that  if  for  some  T  £  R+,  and  ui  and 
u2  such  that  ||PyUi||2  <  oo  and  ||PTu2||2  <  °°>  then  1 1 Pr 1 1 2  <  oo  and  ||PTe2||2  <  in  addition, 
if  Ui,u2  €  /^(R-1-)?  then  £1,62  £  p2(R+)- 

For  this  system,  the  following  assumption  is  made. 

Assumption  7.1  (i)  The  admissible  uncertainty  set  is 

BA  :=  {block-diag{ Ai,  ■  ■  •,  AN}  :  A;  is  causal  and  has  C2-gain  <  1,  for  each  i }; 

(ii)  The  nominal  system  is  time-invariant  and  causal;  in  addition,  given  any  e  >  0,  there  is  a 
function  <f>  :  £|(R+)  X  R+— »R+  such  that  for  all  u  £  £|(R+)  ani d  t  £  R+;  0  =  f’iPtU,  t)  >  t  and 

||(/  -  PHu,t))Gu  -  (/  -  PHn,t))G(I  -  Pt)u\\2  <  e  \\Ptu\\2 

where  Pt  is  the  truncation  operator. 

Defining  the  scaling  matrix  set  V  as  before,  the  following  result  is  essentially  from  [27,  28,  20,  26]. 

Proposition  7.2  (Scaled  Small  Gain  Theorem)  Under  assumption  7.1  (i),  the  uncertain  sys¬ 
tem  depicted  above  is  robustly  BIBO  stable  if  there  is  D  £  V  such  that  the  scaled  nominal  system 
DGD~X  has  C2-gain  <  1;  i.e., 

jf  T 

[  ||  DGD~luf  dt  <  f  \\u\fdt  (82) 

Jo  Jo 

for  all  T  £  R+,  and  u  £  £2[0,P]-  Moreover,  under  additional  assumption  7.1  (ii),  the  condition 
(82)  is  also  necessary  for  robust  BIBO  stability. 
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The  nonlinear  time-invariant  systems  satisfying  assumption  7.1  (ii)  are  said  to  have  fading 
memory  (in  C2  sense)  [30].  It  is  noted  that  the  scaling  treatment  for  the  sufficient  conditions  of 
Theorem  3.5  is  essentially  motivated  by  combining  the  scaled  small  gain  condition  in  this  proposition 
with  the  one  on  £2-gain  characterizations  in  section  2.  However,  the  necessity  of  the  conditions 
in  Theorem  3.5  is  not  as  immediate,  one  of  the  reasons  is  that  the  systems  may  not  have  fading 
memory.  Nevertheless,  only  a  weaker  notion  of  the  fading  memory  about  the  nominal  system  is 
needed  for  the  above  proposition  to  hold  [30],  and  a  large  class  of  the  nominal  plants  considered  in 
this  paper  have  weak  fading  memory  (this  issue  is  discussed  in  a  forthcoming  article). 
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